Size and Location Control of Si Nanocrystals at Ion Beam Synthesis in Thin Si0 2 Films 
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Binary collision simulations of high-fluence 1 keV Si + ion implantation into 8 nm thick Si0 2 films on (OOl)Si 
were combined with kinetic Monte Carlo simulations of Si nanocrystal (NC) formation by phase separation 
during annealing. For nonvolatile memory applications, these simulations help to control size and location of 
NCs. For low concentrations of implanted Si, NCs form via nucleation, growth and Ostwald ripening, whereas 
for high concentrations Si separates by spinodal decomposition. In both regimes, NCs form above a thin NC free 
oxide layer at the Si0 2 /Si interface. This, self-adjusted layer has just a thickness appropriate for NC charging 
by direct electron tunneling. Only in the nucleation regime the width of the tunneling oxide and the mean 
NC diameter remain constant during a long annealing period. This behavior originates from the competition 
of Ostwald ripening and Si loss to the Si/Si0 2 interface. The process simulations predict that, for nonvolatile 
memories, the technological demands on NC synthesis are fulfilled best in the nucleation regime. 



Recently, nonvolatile memory concepts based on nanocrys- 
tals (NCs) embedded in the gate oxide of MOS transistors 
have attracted much interest (]]]. For that aim, NCs have been 
synthesized by a variety of techniques like chemical vapor 
deposition |2], ion implantation |3, 4], and Si aerosol depo- 
sition 0]. Ion implantation followed by thermally activated 
precipitation of the implanted impurity atoms is most com- 
patible with current silicon technology. By low-energy Si + 
ion implantation into thin Si0 2 layers on (OOl)Si, NCs of Si 
were formed a few nanometers above the Si/SiO z interface 
This allows charging of the NCs by direct electron tunnel- 
ing, which is a prerequisite for high endurance and low opera- 
tion voltages 1 6 ] . Further optimization of location and size of 
ion beam synthesized NCs for memory application requires a 
deeper understanding of the mechanisms involved, which de- 
termine (i) the built-up of Si supersaturation by high-fluence 
ion implantation and (ii) NC formation by phase separation. 
For that aim, process simulations were divided into two steps. 
The Si implantation was studied using the binary collision 
code TRIDYN |7], which includes dynamic target changes. 
The phase separation of Si from SiO, during subsequent an- 
nealing has been simulated with a kinetic lattice Monte-Carlo 
code, which describes the thermally activated processes. 

The TRIDYN depth profiles are shown in Fig.[2for 1 keV 
Si + ion implantation into Si0 2 . TRIDYN takes into account 
dynamic target changes due to ion deposition, ion erosion and 
ion beam mixing. The input parameters required by the sim- 
ulation include the displacement and surface binding energies 
of target atoms. The displacement energies E d for both, Si 
and O, were assumed to be 8 eV. This value proved to yield 
satisfactory agreement between earlier TRIDYN simulations 
of ion mixing and experiments 1 8 ] and is also consistent with 
the choice of Sigmund and Gras-Marti 1 9] in their theoretical 
treatment of ion mixing. For the present problem, simula- 
tions with E, varied by factors of 0.5 and 2 did not show any 
significant differences of the Si deposition profiles. The sur- 
face binding energies of Si and O are assumed to vary linearly 
with the surface composition in a way that they balance the 
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Figure 1: TRIDYN depth profiles for 1 keV Si + implantation into 
Si0 2 . For comparison a TRIM profile has been added. 



enthalpies of sublimation and decomposition of Si and Si0 2 , 
respectively I Ki ll IB . 

The Si profile broadens by ion beam mixing, sputtering and 
swelling. For comparison, a Si implantation profile calculated 
by TRIM 1 12] has been added to Fig.^ which is much sharper 
than the corresponding TRIDYN profile. Accordingly, the Si 
peak concentration is significantly overestimated by TRIM. 
Additionally, sputtering has led to a Si enrichment at the target 
surface. 

During subsequent annealing, Si implanted into the thin 
SiO z separates from the oxide phase. In general, the process 
of phase separation is expected to be a sequence of physi- 
cal mechanisms like nucleation, growth, and Ostwald ripen- 
ing of Si precipitates or, at higher Si concentrations, spin- 
odal decomposition and interfacial energy minimization of the 
Si/Si0 2 mixture, respectively. However, these mechanisms 
are the result of a variety of elementary events (like bond 
breaking, diffusional jumps of atoms, chemical reactions etc.) 
that occur in random sequence. Here, they are studied by a 
kinetic 3D lattice Monte Carlo (KMC) method, which is dis- 
cussed elsewhere in detail 1 13]. 

The kinetics of Si atoms is described in a solid host ma- 
trix (SiO,), which is the background or "system's vacuum". 
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Figure 3: Evolution of the width of the denuded zone (a) as well 
as the mean NC diameter and density (b) during annealing for both 
regimes, nucleation and spinodal decomposition. 



Figure 2: Snapshots of KMC simulations (top view and cross- 
section) of phase separation in 8 nm thick Si0 9 on (001) Si dur- 
ing annealing. The simulations start from 1 keV Si + TRIDYN pro- 
files for fluences of (a) 2 x 10 15 cm~ 2 , (b) 5 x 10 15 cm~ 2 , and (c) 
1 x 10 16 cm~ 2 , respectively. Two regimes were identified, "nucle- 
ation and growth" (a) and "spinodal decomposition" (b),(c). Addi- 
tionally, percolation is observed at the highest fluence (c). Atoms 
are colored according to their coordination. The 15 nm scale only 
applies for the lower right corner. 



Thereby, an underlying fee lattice has been assumed, which 
is the most isotropic lattice. Within this host dissolved Si dif- 
fuses and can form precipitates. (The lattice spacing was cho- 
sen such that the correct atomic Si density is obtained.) Ap- 
plying the classical lattice gas model with attractive Si-Si in- 
teraction, the energetics is determined by the nearest-neighbor 
Ising model. The Metropolis algorithm lfl4l1 is used to de- 
scribe the kinetics of the system. 

Si dissolved in the matrix performs diffusional jumps 
from one lattice site to another with the probability 
1/rexp {— E A /k B T}, where E A is the activation energy of dif- 
fusion, 1/t is the attempt frequency and k B T has its usual 
meaning. Due to a time scale normalization, the internal time 
unit of the KMC simulation is a Monte Carlo step (MCS) 
given by Texp{E A /k B T}, thus depending on the tempera- 
ture. Si dimers, trimers and larger agglomerates can form 



and exhibit a binding energy, which is the product of the Si- 
Si bond strength E B by the number of Si-Si bonds. A Si 
atom having n i Si neighbors jumps to an empty neighbor- 
ing site having Si neighbors with a reduced probability 



1/Texpj- E A + ( ni -n f )E B /k B r} 



if rij > n i 



Otherwise, 



the diffusional jump probability remains valid. 

In principle, the bond strength E B for our KMC simulation 
can be determined from the solubility of Si in Si0 2 via the 
detailed balance of Si attachment/detachment at the Si/Si0 2 
interface. (It has to be taken into account that the coordina- 
tion number in the fee lattice is 12 instead of 4 in the case of 
Si lattice.) However, the diffusivity and the solubility of Si 
in Si0 2 are largely unknown. Thus, a direct relation of simu- 
lation and experimental temperatures is difficult, which holds 
for the time scale too. Nevertheless, the path of systems's evo- 
lution towards equilibrium and the regimes of phase separa- 
tion predicted by KMC simulations may improve the process 
understanding substantially. 

In Fig. 12 snapshots of KMC simulations of phase sepa- 
ration are shown for E B /k B T = 2.0. Thereby, the TRIDYN 
profiles of Fig.[2were used as initial Si distributions. The size 
of the simulation cell is 56 x 56 x Aznm 3 , where Az is 8 nm 
plus swelling due to implantation. Besides the image at the 
lower right corner, a quarter of the simulation cell is shown 
only. The simulation cell borders on a fixed (001) layer of 
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the Si substrate to account for the Si/Si0 2 interface. Periodic 
boundary conditions were applied in the oxide plane, whereas 
reflecting boundary conditions were assumed at the surface. 

Two regimes of phase separation are predicted by the KMC 
simulations. A "nucleation and growth" regime is observed 
for 2 x 10 15 Si + cm -2 . Si NCs form by nucleation and grow 
further at the expense of Si supersaturation (Fig. [2] (a), 100 
kMCS). Later on, NCs grow by Ostwald ripening and fi- 
nally dissolve by Si loss to the Si0 2 /Si interface. Above 
2 x 10 15 Si + cm~ 2 a "spinodal decomposition" regime is iden- 
tified, where due to the vanishing nucleation barrier II 511 non- 
spherical, elongated Si structures are formed (Fig. 0(b), 100 
kMCS). At even higher Si concentrations (1 x 10 16 cm~ 2 ), 
above the percolation threshold, the phase separated Si be- 
comes laterally connected (Fig.|2](c), 100 kMCS). This net- 
work of Si does not decay into droplets during longer anneal- 
ing (Fig. |2(c), 4000 kMCS). An electrical charge brought to 
this network can spread over several tens of nanometer, i.e. the 
phase separated Si behaves like a floating gate in a conven- 
tional metal-oxide-silicon transistor. 

Below the percolation threshold at 5 x 10 15 cm~ 2 , the ini- 
tially non-spherical Si structures evolve into spherical NCs, 
which can hardly be distinguished at this late stage from that 
formed by nucleation and growth (see Fig. [2] (a),(b), 4000 
kMCS). The Si0 2 /Si interface acts as an effective sink for 
Si in both regimes, which results in a zone denuded of NCs. 
However, a more detailed consideration reveals some differ- 
ences between the nucleation regime and the spinodal decom- 
position regime. As shown in Fig. [3] a), for the nucleation 
regime (2 x 10 15 cm~ 2 ) the width of the denuded zone is con- 
stant over long annealing times. A similar dependence is ob- 
served for the NC size (Fig.[3](b)). The competition between 
Ostwald ripening and Si loss to the Si/Si0 2 interface keeps the 
mean NC diameter constant over a long period of annealing. 
The Si loss manifests itself in a rapid decrease in NC density 
(Fig. 13(b)). On the other hand, for the spinodal decomposi- 
tion regime(> 5 x 10 15 cm -2 ), the interface minimization of 



the non-spherical Si structures leads to a narrowing of the 
denuded zone (Fig. [3] (a)). Moreover, Ostwald ripening is 
more effective than Si loss to the interface and, hence, up to 
4000 kMCS the mean NC size is increasing as can be seen 
in Fig.[3](b). In contrast to the nucleation regime, the NCs 
become larger and dissolve slower ( Fig. 0(b)). Of course, a 
long lasting annealing results also in this regime in a complete 
dissolution of NCs. 

It should be emphasized that NCs form behind a zone com- 
pletely denuded of NCs (see cross section views of Fig, [2ji. 
Additionally, the distance of the NCs from the interface is 
small enough to allow their charging by direct electron tun- 
neling. This self-alignment of NCs as well as their degrada- 
tion free charging/decharging is crucial for application in non- 
volatile memories |6]. Varying ion implantation energy and 
annealing temperature gives additional control over the width 
of the depleted zone and the NC size, which will be described 
elsewhere flrjll . Using Si implantation profiles predicted by 
TRIM instead of TRIDYN lead to a strong overestimation of 
the percolation behavior during phase separation, which will 
be shown there too. 

In conclusion, two regimes of Si NC formation by 1 keV 
Si + ion implantation into thin Si0 2 and subsequent anneal- 
ing have been found by process simulation. Below a Si flu- 
ence of 2 x 10 15 cm~ 2 , NCs form by nucleation and growth, 
while at higher fluences spinodal decomposition occurs. At 
1 x 10 16 cm~ 2 , percolation leads to a spatially connected 2D 
pattern of Si in Si0 2 . The present KMC simulations pre- 
dict that NC fabrication for nonvolatile memory applications 
should be performed preferably in the nucleation regime. In 
this regime, the width of the denuded zone does not depend 
critically on the annealing time and/or temperature. The syn- 
thesized NCs are small (2.. 3 nm diameter) and of high area 
density (> 10 12 cm -2 ). 
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